Dark halos built of scalar gravitons: numerical study 
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Abstract 



In a previous article due to one of the present authors (YFP), an extension to General 
Relativity, violating general covariance to the residual unimodular one, was proposed. As a 
manifestation of such a violation, there appears the (massive) scalar graviton in addition to 
the massless tensor one. The former was proposed as a candidate on the dark matter in the 
Universe. In a subsequent article (Yu. F. Pirogov, MPLA 24, 3239, 2009; |arXiv:0909.3311| 
[gr-qc]), an application of the extension was developed. Particularly, a regular solution to 
the static spherically symmetric equations in empty space was studied by means of ana- 
lytical methods. This solution was proposed as a prototype model for the galaxy soft-core 
dark halos, with the coherent scalar-graviton field as dark matter. The present report is 
a supplement to the aforementioned article. The statements of the latter are verified and 
visualized by means of numerical analysis and symbolic calculations. The nice validity of 
analytical results is found. 

1 Introduction 

In Ref. PQ, an extension to General Relativity was proposed. The extension possesses the 
residual unimodular covariance, and in line with the massless tensor graviton describes the 
(massive) scalar one. The latter was proposed as a candidate on dark matter in the Universe. 
The theory was further developed in a series of subsequent articles. In particular, in [2] a 
regular solution to the static spherically symmetric equations of extended gravity in empty 
space was qualitatively studied by means of analytical methods. This solution was proposed as 
a prototype model for the galaxy soft-core dark halos, with the coherent scalar-graviton field 
as dark matter. For details, we refer the reader to [U [2]. The present report is a supplement to 
Ref. [2] . The statements of the latter are verified and visualized by means of numerical analysis 
and symbolic calculations. The nice consistency of the qualitative analytical study is found. As 
a by-product, it is found highly plausible that the power series for the regular solution has just 
a finite radius of convergence. 

2 Extended gravity equations 

Here, we shortly remind the results of [2j concerning the static spherically symmetric equations 
of the unimodular extended gravity in empty space. The line element in such a case looks in 
the polar coordinates (r, r, 9, ip) most generally like: 




(1) 
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where a, b and c are some metric potentials. In the static limit, all the variables depend only 
on the radius r, which is still defined ambiguously. To eliminate such an ambiguity impose the 
coordinate condition ab = 1. Defining now the new variables 

A = o = l/6, C = r 2 c (2) 

and supplementing them by the (dimensionless) scalar-graviton field X we get the system of 
the nonlinear differential equations as follows: 
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Here, a dimensionless parameter = 2k 2 /k 2 refers to the Lagrangian of extended gravity, with 
Kg being the mass scale for the ordinary gravity (of the order of Planck mass) and Kh being 
an additional mass scale characteristic of the unimodular extended gravity. On the contrary, 
the parameter > 0, with the dimension of length, arises as a free integration constant. 
Physically, it characterizes an internal length scale for a vacuum solution. 

The first equation above reflects continuity condition in empty space, with the other equa- 
tions being a combination of the continuity condition and three gravity equations. Out of four 
equations, only three are independent and can be chosen at will as a primary system. The 
remaining equation should, generally, serve as a consistency condition. 

At a finite Rh, it is appropriate to choose the scaled distance £ = r/Rh as an independent 
variable. Introducing then t = £ 2 and redefining C as C = tc we get equivalently: 

2 4t{ AC Jt X ) +AC Jt X = 3Cexp(-X), 

2t Jt{ C Jt A ) +C Jt A = & iCexp(-X), 
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This system is the main concern of the present investigation. 



3 Analytical study 

Here, we recapitulate some results of [2] concerning the solution regular at t = 0. Looking for 
such a solution as a power series of t one gets: 
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This representation is, generally, valid at the arbitrary < 1 and may formally be continued 
up to the arbitrary powers of t. But being restricted just to a region of t (see later) it does not 
give apprehension of solution whole. 

To proceed, we restrict ourselves to the case ej, C 1 which is preferred from astronomical 
observations. Under this assumption, decompose formally an arbitrary solution as a power 



series of e\: 



X(t) = f2d n Xn(t), a(t) = J24 n an(t), c(t) = J24 n Cn(t), (6) 

n=0 n=0 n=0 

with the conditions oo = cq = 1. Substituting this decomposition into Q we simplify the latter 
in the leading e^-order as follows: 

d 2 d 2 

dt X * + -^ X » = ex P(-*°)< 
d 2d 2 , v , 

dt ai + 3 t dT 2ai = ex P(- x °)< 

3d ( d ' 

tdt ci+2 df 2Ci = "U ° 

ai - t dt ai - 2t df 2ai = Cl + 5t dt Cl + 2t df 2Cl - (7) 

Clearly, a\ = Xq under the proper boundary conditions. Thus, solving the system of the coupled 
differential equations Q reduces to solving the ordinary differential equation for Xq and c\, 
with the last equation of system serving as a constraint. 
To this end, introducing the new variables 

Z = X - a, a = In (St) (8) 

transform the equation for Xq to the autonomous (not containing explicitly the independent 
variable) form 

^- z Al z= \^- z) - r >- (9) 

Putting then Z = dZ/da reduce the second-order equation for Z to system of two first-order 
equations for Z and Z as follows: 

da 

±i = -\z + \(^ ( -z)-i). (10) 

Such a system is known to be fully characterized by its phase plane (Z, Z). Of the particular 
importance are the exceptional points defined by dZja = dZ J da = 0. In the case at hand, 
there is just one point of this kind, Z = Z = 0, and it belongs to the stable focus type. All the 
trajectories (Z(a), Z(a)) winds around this point approaching it at a — > oo. Moreover, among 

the trajectories there is a unique one (Z(a), Z(a)), with Z remaining finite at a — > — oo. 

In the original terms, this signifies two important properties of the solutions A"o(i). First, 
all of them ripple around the exceptional solution 

X = ln(3i) (11) 

approaching the latter at t — > oo. Second, the regular at the origin solution Xq is unique and, 
supplemented by condition Xq(0) = 0, should look at t — > as given by ([5J with eh = 0: 

Xq = t - -t 2 + -t 3 . (12) 
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This gives the qualitative picture of the looked-for regular solution Xq as a whole. The same 
concerns the regular 5i = Xq. As for the regular ci, it have to approach asymptotically the 
exceptional solution 

ci = 2-ln(3i) (13) 



and be approximated at < t < 1 by the power series 



ci 



-t 2 + -t 3 
10 35 



(14) 



4 Numerical analysis 

Here, we numerically verify and visualize the above analytical statements of [2]. The phase 



plane of system ( 10 ) is shown below: 
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The arrows present the direction field given by dZ/dZ = (dZ / da) / (dZ / da) . The trajectories are 
the lines Z(Z) tangential to the direction field in every point. The bulk of trajectories (dashed 
black lines) possesses Z — > — oo at a — > — oo resulting in the singular at t — > solutions. There 
is just one trajectory, with the finite Z {Z — > — 1 at a — > —oo ), indicated by the solid red line. 
It may be associated with solution Xq regular at t = 0. The exceptional trajectory (Z = Z = 0) 
corresponds to exceptional solution Xq. The picture above explicitly supports the statements 
made in [2]. 

To find the exact form of regular solution we integrate the first equation of system ([7]) 
numerically, with boundary condition taken from (12). This gives 
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Numerical solution for X„ 
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It is well seen that the power series approximates numerical solution very good in interval 
< t < 1. Asymptotically, the regular solution Xq vs. the exceptional one Xq looks like: 
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Clearly, Xq ripples around Xq approaching the latter at t 
approach is though extremely slow. 



oo. The picture shows that this 
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The same concerns c\. Integrating the third equation of system Q numerically we get 
solution as follows: 
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Numerical solution for c x 



At large t, the solution also ripples around the exceptional one approaching the latter asymp- 
totically: 
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Thus, numerical study totally confirms all the qualitative statements about regular solution 
made in [2]. 

At last, to check the accuracy of numerical calculations we input the found numerical solu- 
tions into consistency condition given by the last equation of Q. The absolute (the difference 
of L.H.S. and R.H.S.) and relative (the ratio of L.H.S. and R.H.S. minus unity) errors of calcu- 
lations are shown, respectively, by the blue and red lines below: 

0.003 




-0.002 



-0.003 
0.00015 



0.00010 




-0.00015 



Evidently, the achieved accuracy of numerical calculations is quite nice. 

5 Symbolic calculations 

Here, we examine the exact regular solution to exact system Q in a vicinity of t = under an 
arbitrary Eh- To this end, decompose the looked- for exact solution X(t) as the power series of 
t as follows: 

oo 

x = Y,t n Mth), (15) 

n=0 

with at being some parameters (ao = 0). We proceed similarly with a(t) and c(t)). Implement- 
ing a system of symbolic calculations and following the cyclic perturbative procedure proposed 
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in [2] we get in a reasonable time the 24 terms of decomposition, the first five of them being 
shown below: 



Q'l 
Q'2 

a 3 
0:4 
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(3 + 54), 
(12 + 414 + 35ei) , 

(305 + 15734 + 27354 + 15754) , 
(3774 + 259694 + 681204 + 796754 + 346504) 
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(16) 



(and similarly for a\ and c\). The first three terms above reproduce those given by ([5]). 

Putting now en = we can compare Xq found previously with the present solution, two 
approximations to which being shown below: 
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Numerical solution for X n 



It is seen that the power-series solution X\ £h= o perfectly matches the numerical one Xq up to 
t < 1.7. With n increasing, the approximation clearly improves within this region but worsens 
beyond it. The further analysis reveals the same picture with eh 7^ 0. Fixing a numerical value 
for £/j allows one to calculate much more coefficients a n in a reasonable time. Thus, calculations 
with Sfi = up to n ~ 3- 10 7 make it highly plausible that the area of convergence of power series 
is limited in this case to t < 1.8. It seems that similar statement survives with an arbitrary 
e h <0{\). 



6 Dark halos 

Here, we examine validity of analytical results for rotation velocity and the ensuing dark matter 
profile found in [2] . The circular rotation velocity of a test particle in the spherically symmetric 
metric Q is, generally, as follows: 



v 2 {r) = —a(r) /—In (r 2 c(r)) . (17) 
dr I dr v 



With account for a± = Xq, the respective velocity squared profile (in terms of £ = r/Rh) in the 
leading e/j-order looks like 

vUO/4 = U h (0 = (£/2)dX„(0/de- (18) 
The regular solution Xq results then in: 

u h (0 = I ^ ^ +35^ u s t < 1, (19) 

This analytical approximation vs. numerical result is shown below: 
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Numerical solution for U h (<fj) 




with the numerical Uh approaching unity asymptotically like 
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The regular solution implies the dark matter profile as follows: 

Ph(0/Ph(0) = P h (0 = exp (-Xo(O) , 



where Ph(0) = Qs^k 2 /R\. Analytically, this results in the soft-core halo profile [2]: 



i-e 2 + fe 4 , o < e < i, 
\c\ e » i, 



(20) 



(21) 



with the finite central density Ph(0)- The analytical approximation vs. numerical result is as 
follows: 
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Clearly, the accuracy of analytical study is quite reasonable. 

7 Conclusion 

Summarizing, numerical analysis of regular solution to the static spherically symmetric equa- 
tions of the unimodular extended gravity in empty space totally confirms and somewhat refines 
analytical statements made in pj. As a by-product, symbolic calculations make it highly plau- 
sible that the power-series decomposition of solution is valid just interior to a finite convergence 
radius. We are going to expand the conducted study on the general solutions to the aforemen- 
tioned equations, with a view to refine the application of theory to the galaxy dark halos started 
in 0. 
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